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Vibrational  relaxation  of  OH(X ^ri/,  v=:^-3)  by  O2 

James  A.  Dodd,  Steven  J.  Upson,  and  William  A.  M.  Blumberg 

Geophysics  Laboratory  (AFSO/OPI,  Hanscom  Air  Force  Base.  Massachusetts  01 731 

(Received  8  August  1989;  accepted  3  October  1989) 

Rate  constants  for  OH  (A"  ^H,,  u  =  1-3)  vibrational  relaxation  induced  by  nonreactive  collision 
with  Oi  have  been  measured.  OH(u)  is  created  by  the  H  +  Oj  -•OH(y<9)  +  reaction  in 
an  electron-irradiated  O,,  Hj,  Ar  mixture.  OH(i;)  fundamental  and  first  overtone  IR  emission 
is  observed  using  time-resolved  Fourier  spectroscopy.  Spectral  fitting  followed  by  kinetic 
fitting  of  the  resultant  populations  using  a  single-quantum  relaxation  model  yields  rate 
constants  of  A-,.  _  ,  =  ( 1.3  +  0.4)  X  10~  '  \  A,,  _  2  =  (2.1  +  0;3)  X  10  ‘  A,.„  3  =  (2.9  +  0.8) 

X  10  ^  ' '  (all  units  are  in  cm  Vs).  Our  measu'-ements  are  consistent  with  and  extend  published 
results  on  the  same  system,  as  well  as  predictions  made  by  Schwartz-Slawsky-Herzfeld  theory. 


EXPERIMENT 


INTRODUCTION 

Hydroxyl  radical  (OH )  nlav.s  a  preeminent  rclr  in  both 
the  chemistry  and  photophysics  of  the  earth's  atmosphere. 
Stratospheric  OH  takes  part  in  a  number  of  key  reactions, 
including  those  which  determine  the  concentration  of  the 
halogen-containing  species  involved  in  the  catalytic  destruc¬ 
tion  of  ozone.’  OH  also  mediates  the  relative  densities  of  the 
odd  nitrogen  oxides  ( NO  and  NO, )  and  nitric  acid.  Because 
the  presence  of  OH  is  a  signature  of  numerous  chemical  and 
photochemical  processes,  its  altitude  distribution  is  of  fun¬ 
damental  interest,  and  substantial  effort  has  been^made  to 
develop  better  remote  sensing  techniques  for  OH,’ 

MciiieP  first  identified  the  intense  nighttime  emission  in 
the  near  IR  as  arising  from  OH  overtone  transitions.  The 
reaction  between  hydrogen  atoms  and  ozone  produces  OH 
(A” ’ll,.!  )  in  highly  excited  vibrational  levels,  with  most  of 
the  nascent  product  distribution  in  i’ =  7-9,  i.e.,  21  5(K)- 
26  200  cm  '  above  the  ground  state  (OHt  will  denote  this 
species  with  1:  1).  Unusually  large  probabilities  for  OHt 
multiquantuin  transitions  give  rise  to  the  observed  wide¬ 
band  emission.  Theoretical  and  recent  laboratory  work  has 
addressed  the  dearth  of  accurately  known  rovibralional  Ein¬ 
stein  coefficients  for  OH+  IR  fluorescence,  which  are  neces¬ 
sary  to  determine  populations  frym  observed  emission  inten¬ 
sities. 

^  The  steady-state  vibrational  populations  of  hydroxyl  in 
the  upper  atmosphere  are  determined  by  the  rales  of  various 
processes,  including  the  chemical  reaction^hat  form  and 
destroy  OHt.'IR  fluorescence  from  OHf.  and  collisioiial 
relaxation  of  OHt  by  bath  gas  molecules.  vAt  lower  altitudes 
collisioiial  relaxation  becomes  important,  and  accurate 
modeling  of  the  atmosphere  requires  knowledge  of  the  rel¬ 
evant  rale  constants,  especially  with  the  dominant  molecular 
constituents.  Much  less  is  known,  however,  about  the  vibra¬ 
tional  relaxation  rales  of  OH+  than  about  analogous  rates  of 
stable,  closed-shell  diatomic  molecules,.'Such  as  the  hydro¬ 
gen  halides.  In  this  report  we  addre^  the  relaxation  of 
OH(i'  1-3)  upon  nonreacti’  - ..  ••••tl,  .i.\yg,,ii  ,noic- 

v-ules. 


All  experiments  were  performed  using  the  room-tem¬ 
perature  LABCEDE  apparatus,  which  has  been  described  in 
detail. A  pulsed  35  kV  electron  beam  excites  a  mixture  of 
gases  in  a  5.4  liter,  stainless-steel  vacuum  chamber  ("target 
chamber"),  pumped  by  a  mechanical  forepump.  Differen¬ 
tial  pumping  of  the  filament  source  region  allows  total  pres¬ 
sures  of  up  to  1  atm  in  the  target  chamber.  Calibrated  flow- 
controllers  (MKS  1258A;  1259B)  and  a  0-100  Torr  range 
pressure  transducer  (MKS  Baratron  220BA)  were  used  to 
regulate  gas  flows  and  partial  pressuics. 

OHt  was  produced  by  the  reaction'’ 

H -yO.-OH(t%-9) -hO,.  2.8(  -  11)  cmVs.  (1) 

(Hereafter  rate  constants  are  written  in  abbreviated  form; 
eg,  in  the  above  expression  2.8(  —  11)  denotes  2,8 
X  10  ”,  ]  Ozone  was  formed  using  a  commercial  ozonator 
and  stored  on  a  column  of  silica  gel  at  —  78  °C.  Prior  to  data- 
taking  the  column  was  purged  of  O,  by  pumping.  Argon  was 
subsequently  flowed  (35  seem)  through  the  column  to  en¬ 
train  some  of  the  adsorbed  O,  and  pass  it  into  the  target 
chamber.  Although  we  had  no  direct  monitor  for  the  ozone 
concentration  in  the  target  chamber,  we  estimate  about  10” 
molecules/cm'onthe  basis  of  theOH(i'  =  9)  formation  raie 
( see  below ) .  The  small  flow  of  argon  maintains  this  concen¬ 
tration  ofO,  in  the  target  chamber  throughout  the  day;  only 
a  small  fraction  of  the  ozone  on  the  column  is  actually  u.sed 
in  the  experiment. 

The  other  precursor,  H, ,  is  the  source  for  the  H  atoms 
which  react  with  O,  molecules  to  form  OH+.  Both  ozone 
(0.03%  of  the  total  number  density)  and  hydrogen  (1%) 
are  diluted  in  an  inert  buffer  gas,  in  this  case  argon.  Thus, 
typical  flow  rates  and  pressures  were  as  follows:  Ar,  4000 
seem,  10  Torr;  H, ,  60  seem,  0. 10Torr;O„  trace.  In  addition, 
up  to  250  seem  (0.50  Torr)  O,  gas  was  added  as  a  relaxer  for 
OHt.  Ultrahigii  purity  (99.999%)  argon,  hydrogen,  and 
oxygen  all  were  used  without  further  purification.  Impurity 
present  in  the  added  relaxer  gas  (oxygen)  will  result  in  at 
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most  a  1  %  effect  on  the  derived  rate  constants,  assuming  gas 
kinetic  relaxation  of  OHt  by  the  impurity. 

The  electron  beam  is  pulsed  on  for  300  /ts,  then  off  for 
19.7  ms  to  give  a  1.5%  gun  duty  cycle;  observed  current  at 
the  target  chamber  was  2-3  mA.  The  beam  enters  the  target 
chamber  through  a  pinhole  orifice  and  is  spread  through 
elastic  and  inelastic  collisions  with  the  buffer  gas.  Some  frac¬ 
tion  of  the  energy  lost  by  the  beam  results  in  molecular  exci¬ 
tation.  Intense  fluorescence  is  seen  in  neat  argon,  both  dur¬ 
ing  the  beam  pulse  and  for  several  ms  following  pulse 
completion.  When  small  amounts  of  ozone  or  oxygen  are 
added,  this  fluorescence  remains  intense  during  the  beam 
pulse,  but  is  fully  quenched  50//s  after  the  beam  is  turned  off. 
When  hydrogen  is  also  present,  OHt  chemiluminescence 
dominates  the  IR  emission  following  beam  pulse  termi¬ 
nation. 

A  slow-scanning,  room-temperature  Michelson  inter¬ 
ferometer  is  used  to  observe  the  resultant  OHt  emission  at 
90°  to  the  electron  beam  axis  through  a  6.35  cm  diam  CaF, 
window.  The  central  fringe  pattern  is  focused  by  a  CaF,  lens 
onto  a  1.5  mni  diam  InSb  detector  cooled  to  77  K.  The  detec¬ 
tor  responds  from  roughly  1.3  to  5.5//,  D  *  peaking  near  5//. 
This  region  overlaps  the  fundamental  and  first  overtone 
emission  spectrum  of  OHt,  and  both  sets  of  bands  are  evi¬ 
dent  in  all  of  our  spectral  data.  Spectra  were  obtained  with 
1 5-20  cm  ‘  resolution,  which  is  the  order  of  the  rotational 
constant  B  in  OH.  OHt  emission  was  sampled  at  15-70  //s 
intervals,  such  that  about  80  samples  were  used  to  character¬ 
ize  the  total  IR  emission  decay  curve.  The  detector  output 
was  successively  passed  through  a  matched  preamplifier 
( 10 X  gain)  and  a  PAR  113  preamplifier  with  1  Hz  low- 
frequency  and  10  kHz  high-frequency  rolloff  settings  (50- 
200 X  gain),  and  then  digitized. 

The  data-taking  algorithm  is  described  in  detail  else¬ 
where.'  Helium-neon  laser  fringes  are  used  as  indicators  of 
the  relative  mirror  position.  About  25  duty  cycles  of  fluores¬ 
cence  decay  information  are  averaged  during  each  He-Ne 
A  /2  fringe;  each  fringe  is  scanned  in  0.5  s  in  real  time.  Thus,  a 
two-dimensional  matrix  of  intensity  as  a  function  of  fluores¬ 
cence  decay  time  and  interferometer  fringe  position  is  built 
up  over  time.  Typically,  this  matrix  is  composed  of  80  time 
samples  and  2000  mirror  positions.  Data  are  collected  in  real 
time  by  a  Masscomp  5300  computer  using  16  bit  resolution, 
held  in  memory  until  completion  of  the  scan,  then  trans¬ 
ferred  to  an  Apollo  DN3000  computer.  The  Apollo  is  used 
to  perform  16  K  Fourier  transforms  on  the  raw,  double¬ 
sided  interferograms.  Numerical  phase  correction  methods’ 
are  used  to  generate  real,  time-ordered  fluorescence  spectra. 
Spectral  fits  were  performed  using  a  V  \X  8650  computer 
(see  next  section). 

Accurate  correction  of  spectra  for  interferometer-de¬ 
tector  spectral  response  was  crucial,  since  the  emission  spans 
a  wave  number  range  in  which  the  response  varies  by  more 
than  an  order  of  magnitude.  Calibration  was  accomplished 
using  1000  °C  b!ac^hody  radiation,  which  was  expaiiucu 
with  a  CaFj  telescope  focused  at  infinity  in  order  to  fill  the 
field  of  view  directly  in  front  of  the  interferometer  opening. 
The  calibration  was  performed  under  dry  nitrogen  to  elimi¬ 
nate  the  effects  of  strong  atmospheric  absorbers  ( Hj  O  and 


CO,).  The  curve  obtained  was  insensitive  to  the  precise 
placement  of  either  the  blackbody  or  the  telescope.  Identical 
calibration  curves  were  obtained  on  two  different  days  near 
the  beginning  and  end  of  the  study.  Data  were  corrected 
directly  prior  to  spectral  fitting.  A  corrected  spectrum  is 
shown  in  Fig.  1. 

SPECTRAL  FITTING  OF  OHt  FLUORESCENCE 
SPECTRA 

Observed  spectra  show  clear  evidence  of  OHt  emission 
from  V  =  1-9.  It  is  well  known  from  atmospheric  and  experi¬ 
mental  observations  that  OHt  has  anomalously  large  first 
overtone  transition  probabilities."  Thus,  all  eight  v^i’  —  2 
transitions  were  observed  in  the  overtone  region,  extending 
roughly  from  4200  to  7200  cm  '.  The  fundamental  bands 
with  c'  >  6  are  embedded  in  the  noise  of  our  spectra,  however, 
even  though  they  occur  in  a  region  where  the  detector  is 
most  sensitive  ( 1800  to  2800  cm  ' ).  All  spectra  were  fit  in 
the  wave  number  range  2500  to  7300  cm  ^  '. 

Accurate  spectral  fitting  of  OHt  vibrational  emission 
spectra  requires  knowledge  of  both  rovibrational  frequencies 
and  Einstein  coeflicients.  Spin-rotational  teun  values  for 
t’<IO  accurate  to  a  few  hundredths  of  a  cm  '  have  been 
published."  '"  Mies"  and  Langhoff,  Werner,  and  Rosmus'' 
have  reported  detailed  computations  of  theoretical  transi¬ 
tion  probabilities  for  OHt  rovibrational  lines  for  v'  =  1-9, 
Au  =  1-5.  Very  recently,  Turnbull  and  Lowe"  published 
values  based  on  an  empirical  dipole  moment  function  for 
OH,  while  Nelson  et  al.  reported  values  obtained  through  a 
direct  measurement  of  OHt  transition  probabilities  using 
both  absorption  and  emission  spectroscopy.  At  this  time, 
however,  certain  important  aspects  of  OHt  emission  remain 
controversial,  notably  the  relative  vibrational  band  intensi- 
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Wovenumber  (cm  *) 

FIG.  1  Typical  observed  OHt  emission  spectrum,  corrected  for  spectral 
response.  Note  the  prominent  ^branches  in  the  overtone  rep’^'n  for : '  — 
Kesiduais  lutin  the  ,sj>eciidl  tii  aic  shown  ousel  beneath  ihe  spectrum: 
zero  lines  are  included  for  both  plots.  Statistical  weighting  of  data  points  in 
the  fit  causes  the  larger  residuals  evident  at  larger  wave  number  (see  text ). 
Pressures  used;  10. 1  Torr  Ar.  0. 10  Torr  H  , .  0.00.^  Torr  O, .  c-bcam  param¬ 
eters:  35.0  kV.  2  55  mA  al  the  target  chamber.  0.3  ms  duration,  19,7  ms 
between  pulses.  Spectrum  was  taken  0.3  ms  after  beam  termination. 
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ties.'^  As  discussed  below,  this  issue  does  not  have  a  large 
impact  on  our  analysis. 

In  order  to  spectrally  fit  observed  data,  tabulated  inten¬ 
sities  must  be  multiplied  by  initial-state  populations,  and  the 
lines  then  convolved  with  the  interferometer  line-shape 
function.  Under  our  experimental  conditions  collisions  oc¬ 
cur  on  a  time  scale  (0.1  /rs)  much  shorter  than  the  fluores¬ 
cence  decay  times,  suggesting  a  steady-state  Boltzmann  dis¬ 
tribution  of  rotational  states  at  each  u  level.  The  line-shape 
function  is  conveniently  written  as  a  function  of  two  param¬ 
eters:  the  FWHM  linewidth  A,  and  a  parameter  a  which 
describes  the  line  shape.  The  interferometer  output  over 
2000  fringes  using  an  extended,  monochromatic  helium 
lamp  source  can  be  approximated  as  a  linear  combination  of 
an  unapodized  and  a  triangularly  apodized  interferogram. 
Fourier  transformation  of  such  an  interferogram  gives  a 
spectral  line  described  by  a  linear  combination  of  the  sine 
and  sine'  functions. The  parameter  a  specifies  the  fraction 
of  the  sine  function  in  the  superposition. 

Vibrational  populations  were  determined  through  lin¬ 
ear  least-squares  fits  to  the  spectra,  using  basis  functions  ob¬ 
tained  by  convolving  Boltzmann-weighted  line  strengths 
with  the  instrument  line-shape  function.  The  factors  which 
multiply  each  of  the  nine  vibrational  basis  functions  in  the 
best-fit  result  are  the  relative  populations.  Best  values  for  A. 
a.  and  T aic  obtained  through  nonlinear  fits  to  selected  spec¬ 
tra  for  a  given  day’s  data.  These  parameters  are  then  fixed, 
and  basis  functions  are  computed  for  the  far  less  time-con¬ 
suming  linear  fits  for  the  populations.  A,  a.  and  7’  do  not 
vary  more  than  a  few  percent  over  the  course  of  the  day,  and 
the  variation  does  not  appreciably  affect  the  best-fit  values  of 
the  populations.  Typically,  A  =  16-19  cm  '  (depending  on 
how  many  interferogram  fringes  are  used),  a  =  0.75.  and 
/'=  .100  K.  i.e..  room  temperature.  No  evidence  is  seen  at 
any  time  in  the  fluoreseence  decay  of  a  nonthermal  distribu¬ 
tion  of  rotational  states. 

Spectra  not  corrected  for  detector  response  show  equal 
noise  at  all  wave-number  values;  division  by  the  spectral  cor¬ 
rection  function  results  in  a  noise  level  which  increases  with 
wave  number.  Thus,  points  were  statistically  weighted  as  the 
inverse  of  the  spectral  correction  futiction  in  the  linear  fits, 
placing  more  importance  on  data  near  the  red  end  of  the 
spectrum  where  the  detector  is  most  sensitive. 

Both  the  Mies  and  Langhoff  cl  al.  sets  of  Einstein  coeffi¬ 
cients  yield  good  spectral  fits  of  our  data,  as  measured  by  the 
value  of  resulting  from  the  fitting  program.”  In  addi- 
tioti.  the  derived  rate  constants  did  not  change  significantly 
when  one  or  the  other  set  of  Einstein  coefficients  was  used. 
Because  our  derived  rate  constants  are  insensitive  to  the 
choice  of  Einstein  coefficient  data  set,  the  Mies  and  Langhoff 
cl  al.  data  sets  were  used  interchangeably  in  analyzing  Ollt 
spectral  data. 
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FIG.  2.  Observed  population  decay  curves  forOHtr  —  1-4).  derived  from 
spectral  fits.  Also  shown  is  a  kinetic  fit  to  r  —  1-3  with  v  --  4  as  an  empirical 
feed.  Poorer  fits  result  if  populations  immediately  after  beam  termination 
(time  7-0)  are  included.  At  early  times  rapidly  cascading  upper  r  states 
and/or  H  atoms  produce  conditions  which  violate  the  single-quantum  de¬ 
cay  model  (see  text ). 


smaller  and  have  been  omitted  for  clarity.  Ideally,  this  data  is 
least-squares  fit  to  determine  the  rates  describing  all  OH+ 
single  and  multiquantum  relaxation  processes.  However,  se¬ 
verely  imprecise  and  correlated  kinetic  parameters  result 
when  any  subset  containing  n  curves  is  fit  using  more  than  n 
adjustable  parameters.  Because  we  are  unable  to  determine 
all  of  the  decay  rates,  we  use  the  most  salient  model  which 
incorporates  one  rate  for  each  v  level,  i.e..  that  which  allows 
single-quantum  relaxation  only.  Clearly,  this  model  may  fail 
due  to  two  mechanisms  not  taken  into  account:  ( 1 )  multi- 
quantum  relaxation  and  (2)  chemical  quenching,  which  re¬ 
moves  OHt  from  the  r  ladder  altogether.  We  now  present 
evidence  supporting  the  validity  of  our  simplified  model. 

Little  is  known  about  the  detailed  behavior  of  vibration- 
ally  excited  OHt  upon  collision  with  relaxer  molecules. 
Modeling  studies"*  '’'  of  state-to-state  cross  sections  for 
OH(t’,y)  collisions  with  Ar  atoms,  however,  suggest  that 
multiquantum  relaxation  is  important  for  c  =  9,  but  not  im¬ 
portant  for  V  =  1-4.  These  studies  suggest  that  the  single¬ 
quantum  mechanism  will  model  the  low  v  levels  better  than 
the  high  i:  levels. 

Two  kinds  of  species  which  chemically  quench  OHt 


need  to  be  considered:  those  which  are  flowed  through  the  ^ 
apparatus  continuously,  and  those  which  are  created  by  the 
action  of  the  electron  beam.  In  the  first  category,  only  ozone  q 
reacts  with  OHt  with  a  rate  that  may  be  significant:  q 


OH -t  H,-H -f  H,0,  6.3(-15)cmVs 


KINETIC  MODEL  FOR  FITTING  OF  POPULATION 
DECAYS:  RELEVANT  CHEMISTRY 

The  spectral  fits  produce  a  family  of  population  decay 
curves,  with  =  9.  Figure  2  shows  such  a  family  of  curves 
for  the  four  lowest  levels;  the  populations  in  n  =  5-9  are 


(20  s  '  at  0.10  Torr  H,),  (2) 

OH -(- 0,->  products,  very  slow,  (1) 

OH(n  =  0) -hO,-HO, -1-0,,  6.7(  -  14)  cmVs 

to 

(6  s  'at  0.002  Torr  O,).  (4)  r 

OH(i»l) -tO.-IU-rrChv.  (5a) 
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OH(n>3)  +  0,-0H(t'')  +  O  +  O,,  (5b) 

OH(('  =  9)  +  O,  — products,  2.0(  —  10)  cm  Vs 

(18  000  s  '  at  0.003  Torr  O,). 

(5c) 

Two  different  reports’"  ’'  have  indicated  that  the  reaction  of 
OH(9)  with  O,  is  very  rapid.  The  rate  constants  for  the 
detailed  set  of  OH/O,  reactions  indicated  above  have  not 
been  determined  for  ls,5i'%,8.  However,  new  channels  open 
up  with  increasing  excitation  of  OH.  In  addition  to  the  above 
reactions,  for  instance,  the  channel 

OH(n)  ^  0:-H0,  4- O  (6) 

is  energetically  feasible  for  06.  Since  the  lower  v  levels 
would  be  expected  to  have  longer  reactive  lifetimes,  they 
should  be  more  accurately  fit  by  the  single-quantum  relaxa¬ 
tion  model.  Due  to  its  numerous  reactions  with  OH+,  the 
concentration  of  O,  was  minimized  consistent  with  S/N  re¬ 
quirements. 

Species  which  are  created  by  the  electron  beam  often 
cannot  be  observed  directly,  but  must  be  identified  from  con¬ 
siderations  of  the  reaction  chemistry.  Ions  formed  during  the 
beam  pulse  recombine  with  electrons  on  a  time  scale  which  is 
short  compared  to  that  of  OH  fluorescence  decay.’’  Atoms, 
in  this  case  H  and  O,  pose  a  more  serious  problem.  Hydrogen 
atoms  in  all  likelihood  are  very  efficient  multiquantum  re- 
laxers  of  OH+,  operatitig  through  an  atom  exchange  mecha¬ 
nism.’’  However,  H  atoms  have  a  lifetime  following  beam 
pulse  termination  which  is  limited  by  the  O,  concentration, 

H  -  0,-0Ht  +  O;,  2.8(-ll)cmVs 

(2600  s  '  at  0.003  Torr  O,).  (1') 

Thus,  otic  allows  the  H-atom  population  to  decay  before 
anaU/.ing  the  ()H+  fluorescence  (see  below  ). 

.‘\s  discussed  in  the  literature  OH (9)  reacts  with  O, 
ten  times  faster  than  it  is  produced  by  the  H  -t-  O,  reaction. 
Oroenblatt  and  Wiesenfeld,  "  working  under  pressure  condi¬ 
tions  similar  to  tho.se  reported  here,  found  the  decay  rale  for 
OH  ( r  =  9 )  to  be  independent  of  bath  gas  pressure,  and  iden¬ 
tified  the  rate  associated  with  the  decay  to  OHt  formation, 
not  relaxation.  We  observe  also  that  theOHfi’  =  9)  fluores¬ 
cence  decay  rate  is  independent  of  bath  gas  pressure.  Thus, 
we  identify  the  fluorescence  decay  of  OH(if  =  9)  with  the 
decay  of  H-atom  precursor.  Consistent  with  this  identifica¬ 
tion,  we  obtain  much  better  fits  for  the  kinetic  behavior  of 
OH(t'  =  1-3)  by  first  allowing  OH(i’  =  9)  fluorescence  to 
decay.  In  addition,  one  can  estimate  the  O,  concentration  by 
dividing  the  observed  OH(i'  =  9)  decay  rate  by  the  rate  con¬ 
stant  given  in  reaction  ( 1 ).  The  decay  rate  of  2600  s  '  leads 
to  the  reported  pressure  of  0.003  Torr  O, . 

Unlike  H  atoms,  ground-state  Oi'P)  does  not  react 
quickly  with  any  of  the  molecules  present  with  the  exception 
of  OH: 

O  ■+■  Hj  — OH  -I-  H,  very  slow,  (7) 

O-H  0,-202,  8.0t  —  15)  cm  Vs 

(1  s  '  at  0.003  Torr  O,),  (8) 

0-f0H-02  +  H,  3.9(  -  11)  cmVs.  (9) 


Reaction  (9)  generates  H  atoms  from  ground-state  OH, 
which  in  turn  react  with  O,  uia  reaction  ( 1 )  to  regenerate 
OH  (i’<9).  This  delayed  feed  mechanism  severely  compli¬ 
cates  the  kinetics.  If  high  concentrations  of  O  exist  in  the  cell, 
the  observed  rates  are  lower  than  the  true  relaxation  rates 
due  to  this  delayed  feed.  Anomalously  low  relaxation  rates 
were  observed  when  ( 1 )  the  gun  duly  cycle  was  increased 
from  1.59f  to  3  %,  and  (2)  the  O,  pressure  was  increased 
above  about  0.5  Torr,  at  which  point  the  Stern-Volmer  plot 
of  rate  vs  O.  pressure  takes  on  negaiisc  curvature 

The  effect  of  increasing  the  beam-on  time  is  to  produce 
more  of  all  beam-created  species,  the  steady-state  popula¬ 
tions  of  which  are  complex  functions  of  their  creation, 
quenching,  and  pump-out  rates.  Two  instrunietiidl  time 
scales  of  note  are  the  residence  time  in  the  field  of  \  iew  of  the 
interferometer  (O.I  s)  and  in  the  target  chamber  as  a  whole 
( 1  s).  These  time  scales  are  longer  than  the  OH+  fluores¬ 
cence  decay  (5  ms)  and  duty  cycle  (20  ms)  times.  Thus, 
long-lived  species  may  play  a  significant  role  in  OHt 
quenching.  For  instance,  if  the  OH  concentration  is  too  high 
self-degradation  reactions  become  important: 

OH +  OH(n)-OH +OH(i'),  (10) 

OH +  0H-0-t- H2O,  2.0(  -  12)  cniVs.  (11) 

To  guard  against  quenching  by  beam-created  species,  one 
either  ( 1 )  extrapolates  rates  observed  at  two  or  more  pulse 
widths  to  zero  pulse  width,  or  (2)  demonstrates  that  the 
rates  do  not  depend  on  pulse  width.  The  latter  approach  was 
used  here;  no  dependence  of  the  rates  was  found  as  a  function 
of  pul.se  width  from  250-350 //s.  A  pulse  width  of  300 /ts  w  as 
used  in  most  of  the  experiments. 

The  effect  of  O,  pressure  on  O-atom  density  can  be  un¬ 
derstood  by  considering  the  chemistry  of  the  Ar/H  ./0,/0 . 
system  under  c-beam  excitation.  High  energy  excitation  tif 
this  gas  mixtuie  produces  H  atoms,  as  well  as  O  atoms  in  the 
'/’and  '//states.  H  reacts  with  O,  ria  reaction  ( 1 )  to  pro¬ 
duce  OHt  in  highly  excited  vibrational  states.  Any  0(  '/>) 
present  during  the  e-beam  is  quickly  quenched  after  e-beam 
turn-off.  either  by  reaction  w'ith  H,  to  produce  OHt  in  low  r 
(|)s;4  is  energetically  allowed)  and  high 

0('/») -f  H.-OHt H  3(  -  10)  cmVs 

(960  000  s  'at  0.10  Torr  H.)  (12) 

or  by  collision  with  Ar  atoms’"  or  O,  molecules,'^ 

0( '/))-(- Ar  —  0('/’) -H  Ar,  3.6(  —  13)  cm  Vs 

(115  000s  '  at  10.0  Torr  Ar).  (13) 

0('/)) -f  02-0('P) O2,  2.4(-Il)  cm  Vs 

(370  000)  s  '  at  0.5  Torr  O,).  (14) 

0(  'F)  densities  increase  at  higher  O2  pressure,  due  to  a  fas¬ 
ter  production  rate  from  Oj  and  increased  quenching  ri'a 
reaction  (14).  Thus,  to  avoid  the  delayed  feed  mechanism 
discussed  above,  we  derive  the  OHt/Oj  relaxation  rate  con¬ 
stant  in  the  O2  pressure  regime  where  the  slope  of  the  Stem- 
Volmer  plot  is  linear,  and  the  rates  do  not  depend  on  pulse 
width. 

In  terms  of  the  kinetic  analysis,  clearly  the  levels  which 
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are  fed  directly  by  i^  .ction  (1)  (i>  =  6-9)  cannot  be  de¬ 
scribed  by  our  single-quantum  relaxation  model.  Even  if  a 
single-q"antum  relaxation  mechanism  were  operative  for  all 
i’.  which  .seems  doubtful,  our  model  would  fail  because  the 
observed  fluorescence  decay  for  higher  v  describes  a  forma¬ 
tion  rate;  therefore,  the  fluorescence  for  these  v  levels  cannot 
be  used  to  extract  relaxation  rate  information.  To  analyze 
the  relaxation  kinetics  for  high  i>  a  different  source  must  be 
used. 

Because  of  the  problems  with  the  upper  v  levels,  we  de¬ 
cided  to  concentrate  on  the  lower  v.  In  order  to  provide  a 
(necessary)  feed  term  for  the  highest  v  level  in  the  fit,  we  use 
the  observed  decay  of  a  level  v,  as  an  empirical  feed  term  for 
the  level  directly  oeiow  it.  In  .so  doing  we  remove  c,  from  the 
:.et  of  levels  being  fit.  in  additton  to  all  levels  above  i’,.  In 
order  to  fit  as  matiy  levels  as  possible,  we  choose  Vj  as  the 
highest  level  for  which  the  feed  model  produces  acceptable 
fits.  We  assume  negligible  feed  into  all  v  <  u,  from  all  levels 
above  I'j.  and  that  all  !'<(y  decay  by  the  single-quantum 
mechanism.  We  are  able  to  fit  r  =  1-3  very  well  using  this 
technique,  with  c  =  4  decay  providing  the  empirical  feed 
term. 

The  relevant  kinetic  proces.se.s  included  in  our  model  are 
given  by 

OH(f)  +0,-OH(f-  1)  4  0,,  A(r). 

OH(i')  ^  .V/-OH(c  -  1 )  +  .W, 

OH(c)-0H(f  -  1 )  -  /n.  A 
OH(t’)  -OH(f  -  2)  ^ ,. 

The  rate  constants  k{r)  descrtbe  single-quantum  relaxation 
of  OH((')  by  O,,  the  rate  constants  A  „  describe  relaxation 
by  molecules  other  than  O,  (e.g..  Ar),  and  the  radiative  rates 
A„  are  for  i'  =  i  — n  =  j.  Rates  R  which  do  not  depend  on 
[O,  ]  can  be  grouped  together,  rlz. 

R{\)  ^  A„(l)l.V/l  -V 
/f(2>=--Av,(2)[.V/|  4-T,,. 

R‘{2)  =  R(l)  +  T,„. 

/?(3)  =  A„(3)[,W]  4-/1.,, 

R  '(3)  =  f?(3)  4-/1.,. 

Thus,  the  differential  equations  describing  the  time  depend¬ 
ence  of  the  populations  are  written 

rf  [OH(l)l 
dt 

=  -{/?(!)  -f-A(l)(0,|}(0H(l)l 

4-{R(2)  +A(2)(0,1}(0H(2)I  +/l.,(OH(3)l. 

^  -  {/?'(2)  4- A(2)[0,1}[0H(2)| 
dt 

4-{/?(3)  +A:(3)[02]}[0H(3)], 

-{/?'(3)  f  A(3)lOJ}[OH(3)l 
at 

d[OH(4)) 

dt 

The  A,j  are  about  2  orders  of  magnitude  smaller  than  the 


collisional  relaxation  rates;  the  three-quantum  Einstein  co¬ 
efficients  are  negligible.  The  argon  bath  gas  gives  rise  to  a  y 
intercept  in  the  Stern-Volmer  plot  of  about  one-third  of  the 
maximum  rate  increase  due  to  added  relaxer  O, . 

We  solve  these  coupled  linear  differential  equations  nu¬ 
merically  by  estimating  initial  values  for  the  devay  rates, 
then  integrating  the  equations  over  the  time  of  the  experi¬ 
mental  observations  to  produce  model  curves  for  r  =  1-3.  A 
nonlinear  least-squares-fitting  routine  then  iteratively  ad¬ 
justs  the  rates  to  produce  the  best  least-squares  fit  to  the 
experimental  data,  and  in  so  doing  determines  the  best  val¬ 
ues  of  the  rate  parameters.  Typically  50-80  time  points  for 
each  V  level  are  fit  in  this  manner  for  each  run,  and  several 
runs  are  taken  at  each  relaxer  pressure.  Stern  Vohner  plots 
are  then  constructed  from  the  various  pressures;  a  plot  from 
one  experiment  is  shown  in  Fig.  3.  I  ne  slope  of  the  Stern- 
Volmer  plot  determines  the  rate  constant  for  relaxation  of 
OHt  by  O, .  The  intercept  (typically  1500-2000  s  ')  indi¬ 
cates  the  relaxation  rate  which  is  due  primarily  to  the  buffer 
gas.  Many  days  of  data  were  analyzed,  and  the  rate  constants 
averaged  to  produce  the  rate  constants  reported  herein. 

RESULTS  AND  DISCUSSION 

Rate  constants  derived  for  OHt  relaxation  by  O,  are 
presented  in  Table  1.  We  include  in  the  table  two  previous 
measurements  of  OHt  /O,  relaxation  rate  constants;  both 
are  consistent  with  the  present  measurements.  Jaffer  and 
Smith’"  could  only  obtain  an  upper  bound  for  the  v  =  1  rate 
in  their  flow  tube  experiment.  In  a  more  recent  experiment. 
Rensberger,  Jeffries,  and  Crosley’“  obtained  the  total  rate  for 
relaxation  from  OH (2)  initially  populated  by  a  laser  tuned 
to  various  rotational  lines  in  the  2-0  band.  They  did  not 
determine  the  relative  importance  of  the  single-  vs  double¬ 
quantum  pathways  in  their  experiment.  However,  the  fact 
that  their  value  agrees  with  our  value  arising  from  a  single¬ 
quantum  mechanism  suggests  that  single-quantum  relaxa¬ 
tion  is  more  important  than  double-quantum  relaxation. 

Streit  and  Johnston'"  report  values  for  OH+  relaxation 
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FIG.  3.  Dependence  of  the  decay  rates  for  OH(r  -  1-3)  upon  O,  pressure 
for  one  day's  data.  Nonzero  y  intercepts  are  due  primarily  to  bath-gas-in- 
diiced  relaxation. 
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TABLE  I.  Rate  constants  for  OH  (T  *11,  .n)  +  O, 

-OH(u-  1)  ±0,. 

Rate  constant^  ( 10 

"cm  Vs) 

V 

Prebent  result 

Literature 

1 

1,3  ±0.4 

<4'> 

2 

2.1  -1  0.3 

2.6  +  0.4' 

3 

2.9  ±0.8 

“All  uncertainties  are  Icr. 

''  Reference  28. 

‘  Reference  29.  Measurement  is  the  rate  constant  for  total  collisionally  in¬ 
duced  loss  out  of  i’  =  2. 


by  O,  for  r  =  4-9  in  the  mid- 10  "  cm  Vs  range,  using  a 
kinetic  mode)  which  assumes  single-quantum  relaxation  for 
all  V  levels.  That  work  has  been  strongly  questioned,’'*  both 
because  of  the  complications  described  above  in  the  analysis 
of  the  upper  i\  and  also  because  incorrect  Einstein  coeffi¬ 
cients  were  used.  Finlayson-Pitts  and  Kleindienst ’*  report  a 
rate  constant  for  the  removal  of  OH(9)  by  O,  of  1  X  10  " 
cm  V's,  assuming  the  Mies"  value  for  the /I, ,,  Einstein  coeffi¬ 
cient.  Their  result  is  consistent  with  the  values  reported  here 
for  lower  v.  although  obviously  more  work  needs  to  be  done 
to  determine  rate  con.stants  for  intermediate  v. 

There  have  been  numerous  measurements  of  the  rates  of 
vibrational  relaxation  in  OH  (/I  ’1  ‘  )  by  O;.  *'  For  as 
well  as  for  a  number  of  other  relaxers,  such  rates  are  2  orders 
of  magiiiiudc  larger  than  the  corresponding  values  for 
OH  (A"  Hi, ).  Rensberger,  Jeffries,  and  Crosley’**  suggest  that 
OH(/4)-relaxcr  collision  complexes  have  stronger  binding 
energies  than  their  OH  (A’)  counterparts,  allowing  more  effi¬ 
cient  energy  transfer.  A  more  strongly  bound  complex 
would  be  expected  to  give  rise  to  more  multiquantum  relaxa¬ 
tion  ofOH(i').  German,*-  however,  studying  the  relaxation 
of  OD(/4  "1  ‘  ,1' =  2)  by  O..  obtained  A:(2-*0V 
A(2  — 1)  =  1/10.  a  small  fraction. 

Finally,  we  interpret  our  results  in  the  context  of  the 
Schwartz.  Slawsky,  Herzfeld  (SSH)  theory  "  for  vibrational 
relaxation.  The  large  energy  difference  between  any  of  the 
OH't"  frequencies — v(l)=  3570  cm  i'(9)  —  2237 
cm  ' — and  the  O,  frequency  of  1554  cm  '  makes  resonant 
F-F  transfer  ria 

OHIO  +  OdO)-OH(f'  -  1)  4  0.(11 

infeasible.  The  low  efficiencies  obtained  here  (0.1%)  are 
also  indicative  of  relaxation  which  does  not  use  a  resonant 
F-Fmechanism.  Since  we  are  unable  to  probe  the  vibration¬ 
al  quantum  state  of  the  Oj  product,  we  cannot  distinguish  a 
F-/’  from  a  non-resonant  F-Fmechanism.  For  either  mech¬ 
anism,  however,  SSH  theory  predicts  ( 1 )  rate  constants  for 
single-quantum  relaxation  which  increase  linearly  with  u, 
and  ( 2 )  double-quantum  relaxation  rates  which  are  3  orders 
of  magnitude  smaller  than  the  single-quantum  rates.  Both 
predictions  are  consistent  with  our  experimental  results.  Al¬ 
though  we  obtain  good  fits  to  the  data  using  a  single  quan¬ 
tum  model,  we  cannot  draw  any  conclusions  concerning  the 
precise  importance  of  a  double-quantum  mechanism.  SSH 
theory  fails  when  strong  intermolecular  interactions  occur 


between  the  vibrationally  excited  molecule  and  relaxer  mol¬ 
ecule,  as  is  the  case  at  low  temperatures  (  <  200  K),  or  when 
both  species  are  strongly  polar.  Such  interactions  are  not 
expected  in  the  present  system. 

We  believe  that  our  ability  to  obtain  both  time  and  spec¬ 
trally  resolved  emission  data  of  high  S/N  which  can  be  ana¬ 
lyzed  using  standard  spectral  fitting  techniques  represents  a 
valuable  tool  for  the  deciphering  of  molecular  relaxation 
problems.  In  particular,  our  overall  picture  of  OH'f  decay 
kinetics  complements  the  more  detailed  laser-induced  flu¬ 
orescence  techniques  currently  under  development  for  the 
study  of  high  i;  OHt.  *  We  plan  to  extend  our  studies  of 
OHf  collisional  relaxation  to  other  relaxer  molecules,  in¬ 
cluding  Ar,  NO.  and  N,0. 
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